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We present a theoretical study of the electron-phonon coupling in suspended nanoelectrome- 
chanical systems (NEMS) and investigate the resulting quantum chaotic behavior. The phonons 
are associated with the vibrational modes of a suspended rectangular dielectric plate, with free or 
clamped boundary conditions, whereas the electrons are confined to a large quantum dot (QD) on 
the plate's surface. The deformation potential and piezoelectric interactions are considered. By per- 
forming standard energy-level statistics we demonstrate that the spectral fluctuations exhibit the 
same distributions as those of the Gaussian Orthogonal Ensemble (GOE) or the Gaussian Unitary 
Ensemble (GUE) , therefore evidencing the emergence of quantum chaos. That is verified for a large 
range of material and geometry parameters. In particular, the GUE statistics occurs only in the 
case of a circular QD. It represents an anomalous phenomenon, previously reported for just a small 
number of systems, since the problem is time-reversal invariant. The obtained results are explained 
through a detailed analysis of the Hamiltonian matrix structure. 

PACS numbers: 85.85.+j, 05.45.Mt, 73.21.-b 



I. INTRODUCTION 

The possibility of engineering devices at the nano and 
micro scales has created the conditions for testing fun- 
damental aspects of quantum theoryi, otherwise difficult 
to probe in natural atomic size systems. In particular, 
quantum dots (QD) have largely been considered as a 
physical realization of quantum billiards 2 *^ and meso- 
scopic structures have played an important role in the ex- 
perimental study of quantum chaos^, mainly through the 
investigation of the transport properties of quantum dot& 
and quantum welli structures in the presence of magnetic 
field. However, some extraneous effects can prevent the 
full observation of the quantum chaotic behavior. For in- 
stance, impurities and soft confining potentials may mask 
the chaotic dynamics predicted for some semiconductor 
quantum billiards {e.g., the stadium)^ and the incoher- 
ent influence of the bulk on the electronic dynamics hin- 
ders the observation of the so called eigenstate scars£ in 
quantum corrals^. Furthermore, Random Matrix The- 
ory (RMT) predictions for the Coulomb blockade peaks 
in quantum dots may fail as a result of the coupling with 
the environmenliifl. 

Alternatively, suspended nanostructures are ideal can- 
didates for implementing and investigating coherent phe- 
nomena in semiconductor devices, because, at low tem- 
peratures, they provide excellent isolation for the quan- 
tum system from the bulk of the sample 11,12 . The na- 
noelectromechanical systems (NEMS), in particular, are 
specially suited to study the effects of a phonon bath on 
the electronic states, possibly leading to a chaotic behav- 
ior. Such point is of practical relevance since it bears the 
question of the stability of quantum devicesi^i^, whose 
actual implementation could be prevented by the emer- 
gence of chaos— 

In a recent paper—, we have shown that in fact sus- 
pended nanostructures can display quantum chaotic be- 



havior. In this article we extend such studies and perform 
a detailed analysis of the coupling between the phonons of 
a suspended nanoscopic dielectric plate and the electrons 
of a two-dimensional electron gas (2DEG). The phonons 
are associated with the vibrational modes of a suspended 
rectangular plate (i.e., the phonon cavity) and the 2DEG 
(in the free electron approximation) is confined to a large 
quantum dot (billiard) built on the plate's surface. 

Two different scenarios are considered for the shape 
of the quantum dot: circular and rectangular geometries 
(Fig. [IJ, which yield distinct chaotic features. As for 
the coupling mechanisms, we take into account the de- 
formation and piezoelectric potentials. By performing 
energy-level statistics we show that, for sufficiently strong 
electron-phonon coupling, such electromechanical nanos- 
tructures can exhibit quantum chaos for a large range of 
material and geometry parameters. The resulting spec- 
tral correlation functions, which depend on the geometry 
and location of the center of the QD on the surface of 
the plate as well as on the plate's boundary conditions 
(free or clamped), are those expected from the Gaussian 
Orthogonal Ensemble (GOE) or Gaussian Unitary En- 
semble (GUE) of the Random Matrix Theory (RMT)±1. 
We present a detailed explanation for the occurrence of 
such different statistics distributions. Noteworthy are the 
results for the circular QD, since in this case the GUE 
statistics can be obtained in spite of the fact that the 
system is time-reversal invariant. By investigating the 
influence of material and geometrical parameters on the 
unfolding of chaos, we indicate the conditions for its ex- 
perimental observation. 



II. THE SYSTEM HAMILTONIAN 

The full Hamiltonian of the problem is composed of 
three parts: the phonons, the electrons and the electron- 
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where 



FIG. 1: Schematics of the suspended nanoelectromechanical 
structures depicting the cases of a circular and a rectangular 
quantum dot on the surface of a suspended dielectric plate. 



phonon interactions, which are formulated in the se- 
quence. 
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Likewise for Y n {y). The signs in Eq. © are positive (neg- 
ative) for the FF (CC and CF) boundary conditions 21 . 
The A: m 's are solutions of cos^L^] cosh[/cLj;] = 1. 

Under given boundary conditions, the Rayleigh-Ritz 
method is used to obtain the coefficients A^ nn of Eq. 
(|2*|) and the eigenfrequencies to a corresponding to the 
eigenmode u Q (r) of the plate. It is done by impos- 
ing dU/dA mn = to the energy functional U = 
J dxdy [K,{x, y) — V(x, y)^^, where the kinetic and strain 
energies are written as 



A. Phonons 

At temperatures below 1 Kelvin the acoustic phonon 
mean free path in SiN (silicon nitride), for instance, can 
be as large as 10 /j,mi&. This implies that a plane wave 
acoustic phonon propagating through a suspended meso- 
scopic system whose dimensions are much smaller than 
this mean free path hits the boundaries many times dur- 
ing its expected lifetime, giving rise to standing waves. 
Thus, the phonons in such systems can be described in 
terms of the plate's normal modes of deflection, instead 
of the plane wave phonon description that is more ap- 
propriate for bulk systems. Therefore, in the following 
we associate the phonons with the vibrational modes of 
the suspended mesoscopic system. In addition, at low 
temperatures the semiconductor can be treated as a con- 
tinuum elastic material due to the large wavelength of 
the phonons. 

To obtain the long wavelength vibrational modes of 
the plate we use the Classical Plate Theory (CPT) 
approximation^. The CTP describes adequately the 
vibrations of a plate whose thickness is much smaller 
than its lateral dimensions, which is the characteris- 
tic of our NEMS. The deflections of a plate lying in 
the (x, y) plane are thus described by a vector field 

U(r) i + V(y) ] + W(r) k exp(— iuit) of components 

dW dW 
U(x,y,z) = -z—, V(x,y,z) = -z— , (1) 

(2) 



W(x,y) = A mn X m (x)Y n (y). 



In Eq. (J2J, W(x,y) is written in terms of the one- 
dimensional transverse modes X m and Y n , which are the 
solutions of the Bernoulli-Euler equatio n 19 ; 20 under the 
appropriate boundary conditions. Considering that each 
of the four sides of the plate can be either clamped (C) 
or free (F) (corresponding to the Dirichlet or Neumann 
boundary conditions, respectively), we have 



X m {x) — sin[fc m a;] ± sinh[fc m a;] 

+C{cos[fc m a;] ± cosh[fc m x]} 



fC(x,y) 
V(x,y) 
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Here p 2D denotes the two-dimensional density, v is the 
Poisson constant and D is the rigidity constant. 

In the CPT approximation, the most important motion 
is that in the z direction, given by W(x, y), whereas the 
displacements along the x and y directions are described 
approximately by the first term of an orthogonal basis 
expansion. Therefore, an arbitrary transverse motion can 
be expanded in the basis of the orthonormal vibrational 
modes W a (x,y), for which 



(W a {x,y)\W (x,y)) = S a p, 
J2w a (x,y)W a (x',y')=S(x-x')5(y 



(7) 



The U a and V a components are only approximately or- 
thonormal. As an illustration, we show in Fig. [5] the 
vector field components (U a ,V a ,W a ) of the SA2 eigen- 
mode (i.e., the second symmetric/ antisymmetric eigen- 
mode) for the {FFFF} boundary conditions. Hereafter, 
we refer to the set of boundary conditions of the plate, 
either C or F, as {A1A2Y1I2}, in accordance to Fig. 1. 

An arbitrary vibration field of the cavity is written in 
terms of its deflection modes a as 



u(r,i) = 

a 

x [U a {r)i + V a (r)j+W a (r)k 



(8) 



(3) 



together with the normal coordinates Q a (t) = 
Q a exp[— iu> a t]. In writing Eq. ©, we have taken into 
account that the modes a are real. 

To provide the same level of description for the elastic 
and electronic degrees of freedom of the electromechani- 
cal nanostructure, we perform the canonical quantization 
of the vibration field given by Eq. (jHJ. As a result, we 
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But if we write {X a {t) = Q a (t) + Q\ ) 
W(x,y,t) = J2 X <-(t)W a , 

a 

n z (x,y,t) = pVj2(-^a)X a (t)W a = J2' p c l (t)W a , 

a a 

(10) 

the commutation relation © yields 

[W(x 1 y)^ z (x,y)]=Y}X a . l Vp\W a W p . (11) 

a,/3 

Then, by requiring that \X a ,Vfi\ — ifo5 a ,p one can use 
Eq. Q) to show that Eq. © is satisfied. Thus, X a 
and Va are canonically conjugated operators, satisfying 
[X a ,Xp] = [V a ,V p ] = as well. 

The normal coordinates are now the quantum mechan- 
ical operators Q a (t) and Qa(t), which are used to define 
the dimensionless number operators 



2Vpuj a ^ f 



2Vpui 



(12) 



From the previous commutation relations it can be shown 
that [a a (t),a^{t)\ = 5 a ^ and [a a ,a a ,] = [0^,0^] = 0. 
Therefore, aj, and a a are the creation and annihilation 
operators of the phonon deflection modes a and the vi- 
bration field operator u(r, t) is 



km + aim 

^ ^2 Vp u a /h 



U a (v)i+V a (r)j+W a {r)k 



(13) 



B. Electrons 



FIG. 2: Contour plots for the cavity deflection mode SA2 (the 
second symmetric/antisymmetric eigenmode) for the {FFFF} 
boundary conditions, calculated at the surface of the plate: 
(a) W(x,y), (b) U(x,y), and (c) V(x,y). The amplitudes and 
lateral dimensions are given in arbitrary units. 



associate the classical field u(r, t) with the quantum op- 
erator u(r, t), which must satisfy the equal-time commu- 
tation relation [iij(r, t), TTj(v', £)] = iHS(r — r'), with the 
conjugate momentum operator 7r(r, t) — mdtii. Particu- 
larly, for the k component of the field we have 



We consider the free electron approximation and as- 
sume the electrons to be completely confined to a nar- 
row quantum dot, forming a quasi-2DEG of thickness 
d. The normalized electronic eigenstates are written as 
</> Kj7 (r) = tp K (x, y)y/2/dsm[jTrz/d]. Due to the quasi- 
2D assumption, the electrons always occupy the lowest 
state in the z direction, so in our calculations we set the 
quantum number 7=1. 

For the rectangular QD of sides L x and L y , we have 



Vn(x,y) = 



2 . pirx . . qiry 

sm[— ]sm[— ], (14) 



with k = (p, q) and p, q assuming positive integer values. 
The corresponding eigenenergies are 



E, 



K,7 — 1 



TT 2 h 2 

2m e 



L 2 X L y 



1 



(15) 



[W(x,y),ir z (x,y)] = ihS(x - x)S(y - y). 



(9) where m e is the effective electron mass in the QD. 
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For the circular QD of radius R, we have 
J|i| {u-lvTi) exp[H9] 



y/nR\J\i\ +1 (ai v )\ 



(16) 



with k = (/, v), I — 0, ±1, ±2, . . . and a/„ the v-th. root of 
the Bessel function of order |Z|. Here, the eigenenergies 
are 



we have from Eq. that 



V, 



DP 



\^A a mn \ drz <f>* K „ (XlY n + XrXl) <t>*> ■ 



(21) 
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(17) 



C. Electron-phonon interactions 

The electrons interact with the lattice vibrations 
through different mechanisms, depending on the char- 
acteristics of the solid and the temperature. In addition, 
from the theoretical point of view there can be several ap- 
proaches to describe the coupling between electrons and 
phonons. Next we formulate the electron-phonon inter- 
action terms that are more relevant to our problem. 



1. Deformation potential - DP 

At low temperatures only the long wavelength acoustic 
modes are populated and the semiconductor can be de- 
scribed by the continuum approximation. As a result of 
the cavity deflections, local volume changes take place, 
thus modifying the lattice constant and the electronic 
energy bands. In first order, such volume changes are 
due to longitudinal (compressional) acoustic modes and 
the scattering potential acting on the electrons is propor- 
tional to A(r) = V-u(r). Therefore, the Hamiltonian for 
the DP interaction is 



Hdp — 



C DP I dr #t(r)V • u(r)#(r) 
Jv 



V a k" k' v 

(18) 

with Cpd denoting the deformation potential constant 
for the material. \I>(r) = ~}Z K b K (j) K {r) and ^(r) = 
S K bl-(f)* K (r) are the electron field operators and b K (bj.) is 
the fermionic creation (annihilation) operator satisfying 
usual anti-commutation relations. 

The integral is performed over the volume T> compris- 
ing the 2DEG. In the above equation, V®£, K , is given by 



rf;,';;„, = I drfanV- {u a i+v a j + w a k)<>,,. run 



Since 

V • (UJ + V a j + W a k) = - s (^ + **° ]l (20) 



dx 2 



dy 2 



2. Piezoelectric potential - PZ 

In piezoelectric materials, the acoustic lattice vibra- 
tions produce polarization fields that act back on the vi- 
brational modes. The result is a set of coupled equations 
for the acoustic and polarization fields. However, taking 
into account the difference between the sound and light 
velocities, such equations can be decoupled, yielding the 
following electric field in the semiconductor— 



J - J * y z i & X Z i X y ) • 



(22) 



gX4 and e are, respectively, elements of the piezoelectric 
and dielectric tensors. Expression l|22|l is obtained taking 
into account the cubic symmetry of the lattice. Further- 
more, from the CPT approximation, the strain tensor 
elements e xz = e yz = and 



dU dU 
dy dx 



d 2 W 
dxdy 



(23) 



Therefore, for a given transverse mode a, the resulting 
electric field is perpendicular to the plane of the cavity 



= 2A(z) T 



(24) 



with A(z) = d(2d-z)/2. 

The potential energy of the electrons can be written as 
— e J ~Ei pz ■ d\, leading to 



,egi4 



e A{z) ^iW^J aa+ai] ^xdy-- 



(25) 



Finally, we write down the PZ electron-phonon Hamilto- 
nian as (Cpz = 2egn/e) 



h V pz 



2Vp 



with 



(26) 



PZ 



]T A« nn I dr A(z) ^„X' m Y' n 0„, . (27) 
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D. The full Hamiltonian 

The total Hamiltonian of the system, when both the 
DP and PZ interactions are included, is 

H = H e i + Hph + Hel-ph 

= E R b\b K + ^ (h a + - J Hu a + Hdp + Hpz- 

(28) 

The basis in which H is represented is constructed as 
the product of the one-electron state \<p K ) with the multi- 
phonon state |m, ri2, TI3, njv). Here, n a = 0, 1, . . . , n 
denotes the number of phonon quanta in mode a, with 
maximum population set by n. A total of N distinct 
phonon modes are considered. The values of n and iV 
are chosen to be compatible with the thermodynamics 
of the system. At low temperatures (below IK) n is of 
the order of a few tens, the average phonon occupation 
number. On the other hand, N ranges from ~ 10, at 
the lowest temperatures, up to ~ 30 at the highest ones. 
Hence, in the numerical calculations we set n = 20 and 
N — 21. It has been verified, however, that by varying n 
and N through a considerably wide range does not alter 
our main results. 

A typical basis vector is written, for a given n = 
(ni,n 2 , • . . ,njv), as 

| K ;n)H0 K >®n4^(4r|O}. (29) 

For the diagonalization procedure, we energy-sort the ba- 
sis set up to a maximum energy value. The diagonaliza- 
tion is then performed with such set of vectors. Ob- 
viously, the energy of each basis state is given by the 
sum of the electron and the phonon energies, E Kn — 
E K + J2 a ( n a + l/2)fkj a . E K comes from either Eq. <(T5|) 
or Eq. ijTTI) . depending on the specific geometry of the 
2DEG. For the formation of the original basis set 10 5 
levels were taken into account, however, the diagonaliza- 
tion of H is performed in the truncated basis that varied 
from 3 x 10 3 to 15 x 10 3 basis states. It is important to 
notice that the proportion of different phonon states |n) 
to electron states \<j> K ), comprising the truncated basis, 
ranges from several tens to about a hundred, depend- 
ing on the details of the NEMS. That is, the number 
of phonon states taking part in the calculations is much 
larger. 

III. ENERGY LEVEL STATISTICS 

To determine whether the electron-phonon interaction 
generates chaos in the NEMS, we consider the standard 
approach of looking into the statistical properties of the 
system eigenenergies 5-23 . For completeness, here we give 
a brief summary of the main ideas. A general technical 



overview can be found in Ref. [24| . whereas a very in- 
structive discussion is presented for a particular case in 
Ref. ||. 

Consider the ordered sequence {Ei, E%, . . .} of eigenen- 
ergies of an arbitrary quantum mechanical problem. The 
cumulative spectral function, counting the number of lev- 
els with energy up to E, is written as 

V (E) =J2 @ ( E ~ E ")- ( 3 °) 

n 

In principle, we can always separate rj(E) into smooth 
(average) and oscillatory (fluctuating) parts, so that 

— ^smooth 

(E). (31) 

The smooth part is given by the cumulative mean level 
density 2 ^. 

To make the analysis independent of the particular 
scales of the spectrum, one can use the so called "un- 
folding" procedure^ 3 .. It allows the comparison of the 
results obtained from any specific system with the pre- 
dictions of the RMT— . The unfolding is done basically 
by mapping the sequence {Ex, E2, ■ ■ .} onto the numbers 
{si, s 2 , . . .}, where 

&n — ^/smooth (En). (32) 

In the new variables, the cumulative spectral function 
simply reads fj(s) = s+fj osc (s), so that the smooth part of 
rj has unity derivative. Hence, for our statistical studies 
we consider the resulting sets {si, S2, . . .}. 

In this work we calculate two of the most used spec- 
tral distributions 2 ^: the nearest-neighbor spacing distri- 
bution, P(s), and the spectral rigidity A 3 (^). The P(s) 
distribution probes the short scale fluctuations of the 
spectrum. It corresponds to the probability density of 
two neighboring unfolded levels s n and s„+i being a dis- 
tance s apart. A3 (J) is an example of a distribution that 
quantifies the long scale correlations of the energy spec- 
trum. It measures the deviation of the cumulative num- 
ber of states (within an unfolded energy interval I) from 
a straight line. Formally 

Ag(0 = y (mia {A!B} jf ° ds [ij(s) - As - B] 2 ^J , 

(33) 

where (•) denotes the averaging over different possible 
positions so along the s axes. The parameters A and B 
are chosen to minimize [fj(s) — As — B] 2 in each corre- 
sponding interval. 

The RMT predicts three different classes of Gaus- 
sian ensemblesii^, having distinct P(s) and A3 (2): the 
Gaussian Orthogonal Ensemble (GOE), the Gaussian 
Unitary Ensemble (GUE), and the Gaussian Sympletic 
Ensemble (GSE), constituted by matrices whose elements 
are random and obey certain Gaussian-like distribution 
relationship. Furthermore, these ensembles are invari- 
ant under orthogonal, unitary and sympletic transfor- 
mations, respectively. Bohigas et ali^i conjectured that 
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the spectrum fluctuations of any quantum chaotic sys- 
tem should have the same features of one of such three 
cases. This proposal has been firmly established by theo- 
retical and experimental examinations^Si24. When spin 
is not involved, it is expected that the spectrum statis- 
tics of a chaotic system is similar to that obtained from 
the GOE (GUE) if it is (is not) time reversal invariant 
(TRI). However, there are exceptions to this rule, con- 
sisting of a special class of TRI systems with point group 
irreducible representations, which does exhibit the GUE 
statistical* 2 ^. Until recently^, the only family of sys- 
tems known to show this anomalous behavior was formed 
by billiards having threefold symmetry, implemented ex- 
perimentally in classical microwave cavities22i2ii22. 

For regular (integrable) systems the resulting statistics 
follow Poisson P(s) = exp[— s] and linear As(l) =1/15 
distributions 2 ^. For GOE and GUE, P(s) is described 
with high accuracy by the Wigner distributions 2 ^ 



P(s) = | S exp[-J S 2 ] (GOE), 

32 4 
P(s) = -=s 2 exp[ s 2 ] (GUE). 




Y 



(34) 



FIG. 3: The distinct positions, A, B, C, D, and E (on the 
plate) used as centers for the quantum dot. The dashed 
lines represent the possible symmetry axes for the plate 
phonon modes a, depending on the boundary conditions 
{X 1 X 2 Y 1 Y 2 }. 



Finally, Aa(i) can be approximated by the expressions 
A 3 (/) = i (\n[2nl] + 7 - j - ^ (GOE), 



(GUE). (35) 



Here, 7 = 0.5772 ... is the Euler constant. 

To characterize our nanostructures, we compare the 
numerically calculated distributions P(s) and A 3 (Z) with 
the corresponding analytical expressions for the regular 
and chaotic cases. Very good statistics are obtained using 
2000 up to 2500 energy levels. 



(3 dp and /3pz), the stiffness tensor elements cn,c\2 and 
C44, the mass density of the cavity and the in-plane elec- 
tron effective mass. As for the geometrical parameters, 
we also varied the size and aspect ratio of the dielectric 
plate, the area of the QD, the thicknesses of the plate (<5) 
and of the 2DEG (d). More interestingly, however, we 
considered different positions for the center of the QD 
(shown in Fig. which produces distinct chaotic be- 
haviors. 

The most representative results will be presented 
throughout this section. A detailed analysis is found in 
Section [V] 



IV. RESULTS 



Circular 2DEG 



We have applied the previous analysis to the eigenen- 
ergies of our suspended NEMS, considering a wide range 
of material and geometrical parameters, and it was found 
that chaos emerges in the system for a sufficiently strong 
electron-phonon (el-ph) coupling. Although the phe- 
nomenon proved to be quite robust with respect to vari- 
ations of physical dimensions, boundary conditions and 
basis size, it was observed that the chaotic features de- 
pend on some material parameters, like the electronic 
effective mass and the el-ph coupling constants. The 
materials used to model the NEMS comprise an AlAs 
dielectric phonon cavity and an Alo.5Gao.5As quantum 
dot, where the 2DEG is formed. This choice takes ad- 
vantage of the very small lattice parameter mismatch in 
the interface as well as the large electronic effective mass 
of the X valley in AlGaAsS. 

In our investigation we varied the DP and PZ inter- 
action strengths (by means of the multiplicative factors 



Here we present a detailed analysis for the spectral 
statistics of the NEMS containing a circular quantum 
dot. Unless mentioned otherwise, the system comprises 
a QD of radius R = 450 nm and thickness d — 5/5 on 
the surface of a square phonon cavity of sides L — 1 /im 
and width 6 = 40 nm. 

Once enough electron-phonon coupling is assured, reg- 
ular or chaotic spectral features will emerge depending 
on the interplay between the symmetries of the cavity 
phonon modes and the electronic wavefunctions. In this 
respect the boundary conditions of the phonon cavity 
(i.e., the dielectric plate) and the localization of the cir- 
cular 2DEG play a crucial role. In order to systemati- 
cally investigate this effect we make use of the scheme 
presented in Fig. [3J Slight displacements of the QD 
out of the center of the plate suffice to generate dif- 
ferent spectral features. So, the relative coordinates 
(—0.5 < x,y < 0.5) used in the calculations are: A = 



7 



(0,0), B = (0.05,0.05), C = (0.05,0.025), D =J0.05,0), 
and E = (0,-0.05). Figure H shows P(s) and A 3 (Z) for 
cases A,B,C and D in the {FFFF} phonon cavity, taking 
into account only the DP interaction, with Pop = 10. 
For A, the spectral statistics indicates a regular dynam- 
ics, but in B the occurrence of quantum chaos is clear and 
the level distributions are well described by the predic- 
tions of GOE random matrices. The same occurring for 
D. The more interesting case, however, is C, for which 
the statistics belongs to the GUE class, although the 
system is time- reversal invariant. The same behavior is 
obtained for a phonon cavity with {CCCC} boundary 
conditions^. The reasons for obtaining GUE statistics 
in this time-reversal invariant system will be discussed in 
Section El 
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FIG. 4: Energy-level statistics for the nanostructure with 
only the DP interaction, for Pdp = 10 and four different posi- 
tions (A,B,C, and D) for the center of the QD (refer to Fig. 0. 
The cavity boundary conditions are {FFFF}. The symbols + 
represent the numerically calculated results. The curves indi- 
cate the expected behavior for regular (solid), chaotic GOE- 
type (dashed), and chaotic GUE- type (dot-dashed) systems. 



It is also instructive to look at the evolution of the spec- 
tral statistics as a function of the position of the QD. The 
effect is illustrated by the A 3 (I) statistics in Fig. for 
the {CCCC} plate and the parameters of Fig. 0] Leaving 
A along the S x axis the statistics evolves from regular to 
GOE, at (0.02,0), passing through a mixed behavior at 
the locus (0.005,0). Proceeding perpendicularly to the 
S x axis, the statistics evolves from GOE towards GUE, 
(0.02,0) -> (0.02,0.005) -> (0.02,0.01) = C, going back to 
GOE when the S x - y axis is reached at (0.02,0.02). 

By the examination of several different scenarios we 
were able to classify the general behavior of our system. 
Table I summarizes the results obtained for the center of 
the circular 2DEG located at points A,B,C,D and E with 
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FIG. 5: As(l) x / for different positions of the center of the 
circular QD. The cavity boundary conditions are {CCCC} 
and (3dp = 10. 



either the DP or the PZ interaction taken into account. 
Furthermore, we have considered a comprehensive set of 
boundary conditions, which are representative of all pos- 
sible combinations of the Dirichlet and Neumann condi- 
tions for the phonon cavity, thus, producing distinct sym- 
metry axes for the phonon modes: {CCCC}, {FFFF}, 
{CCFF}, {CFCF}, {CCCF}, and {CFFF}. From Ta- 
ble I, we see that the chaotic behavior is determined by 
the overall (or global) symmetries of the NEMS, that is, 
the one that results from the joint combination of the 
boundary conditions of the phonon cavity and the po- 
sition of the QD. For instance, if the boundary condi- 
tions are {CCCC} and the QD is located at D, the S y , 
S x - V and S x + y are not symmetry axes for the coupled 
electromechanical system. According to Table I, if the 
present NEMS has: (i) four symmetry axes (position A 
for {CCCC} and {FFFF} plates), the statistics indicates 
a regular (integrable) problem; (ii) two symmetry axes 
(position A for {CCFF}), then the statistics corresponds 
to the uncorrelated superposition of two distributions of 
the GOE type^i; (hi) one symmetry axis, the results are 
those of GOE; and finally (iv) no symmetry axes at all, 
the spectrum exhibits the GUE statistics. 

The boundary conditions determine not only the po- 
sition of the center of the 2DEG at which regular, GOE 
or GUE statistics are obtained, but also the intensity of 
the quantum chaos. This happens because the electron- 
phonon coupling depends on the phonon energies, which 
vary according to the boundary conditions. The higher 
the phonon energy, the stronger the electron-phonon in- 
teraction, thus, leading to spectral fluctuations that are 
more faithful to the typical chaotic features. The ener- 
gies of the phonon modes decrease according to the fol- 
lowing sequence: {FFFF}, {CCFF}, {CCCC}, {CFFF}, 
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TABLE I: Symmetry axes and spectral statistics for points A,B,C,D and E, as defined in the text, for a circular 2DEG. Here, 
2 GOE means that the statistics can be described by the uncorrelated superposition of two GOE distributions. 



Boundary Conditions 


Symmetry axes 
for the phonon modes 


A 
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r> 
13 




JJ 


1/ 

Hj 


{FFFF} 


S-jj, Sy, Sx — y, Sx + y 


Regular 


GOE 


GUE 


GOE 


GOE 


{CCCC} 


S x , Sy, Sx — y, Sx + y 


Regular 


GOE 


GUE 


GOE 


GOE 


{CCFF} 


Sx i Sy 


2 GOE 


GUE 


GUE 


GOE 


GOE 


{CFCF} 


Sx — y 


GOE 


GOE 


GUE 


GUE 


GUE 


{CCCF} 


Sy 


GOE 


GUE 


GUE 


GUE 


GOE 


{CFFF} 


Sx 


GOE 


GUE 


GUE 


GOE 


GUE 



{CCCF}, and {CFCF}. The energies for the first five 
boundary conditions are similar, and quantum chaos can 
be observed for essentially the same values of the inter- 
action strength j3. For the {CFCF} case, however, the 
phonon energies can be one order of magnitude smaller 
than the ones for the other cases, requiring larger values 
for the parameters (3dp and (3pz (approximately 3 times 
larger) . 




FIG. 6: The calculated spectral rigidity for various DP cou- 
pling constants: f3op~i (open circle), 3 (filled circle), 5 (x) 
and 10 (+). The system corresponds the {CCCC} plate with 
the QD located at position C. The curves represent the regu- 
lar (solid), GOE (dashed) and GUE (dot-dashed) cases. 

It is, however, important to notice that, regardless of 
the geometry adopted, the nanostructure shows a regular 
spectrum for the bare parameters of the reference materi- 
als. We present in Fig. 0the dependence of the spectral 
rigidity As(l) on the electron-phonon coupling strength. 
For the locus C of the {CCCC} phonon cavity, we take 
into account only the DP interaction, with Pdp = 1,3, 
5, and 10. As (3 increases, the calculated statistics grad- 
ually converges to the GUE prediction. Note that the 
numerical calculations are never well fitted by the GOE 
distribution. The inclusion of more basis states does not 
alter the observed results. At this point we observe that 
the strong el-ph coupling regime (/3 > 3) can be achieved 
by using different materials. For instance, aluminum ni- 
tride (A1N) is a strong piezoelectric semiconductor, with 
£33=1.5 C/m 2 , that is currently been used to produce 
nanomechanical resonators^. The piezoelectric constant 
for GaAs is M =O.16 C/m 2 

Next, we summarize the effects of the geometrical and 



material parameters on the chaotic behavior. Irrespec- 
tive of the boundary conditions, when the QD radius R 
decreases to less than one third of the plate side L, the 
system starts to become regular and, for about L/Rm 5, 
the nanostructure presents no clear signs of chaos in its 
spectrum. This is illustrated in the top panel of Fig. [7J by 
the A 3 (Z) statistics calculated for case A in the {CFFF} 
plate. On the other hand, an important physical param- 
eter for the occurrence of chaos is the in-plane electron 
effective mass m*. It is shown in the bottom panel of 
Fig. [7Jthat chaos arises as m*/m e is increased. Indeed, 
for to* > 0.6 the system is clearly chaotic, becoming reg- 
ular for to* < 0.2 TO e . As before, similar results hold for 
other boundary conditions. A weak dependence on both 
the density and the value of the stiffness tensors cu and 
C44 is also observed. In essence, lighter and softer mate- 
rials favor the appearance of chaos. As for the size of the 
dielectric plate, chaos is favored by short and thin plates. 
The former should be expected because the phonon ener- 
gies increase as the area of the plate decreases. However, 
structures much smaller than the one here considered do 
not present a significantly higher tendency to chaotic be- 
havior. 

So far, we have considered only the DP or PZ inter- 
actions acting individually. When acting together, the 
spectrum statistics of loci B and D change from GOE 
to GUE. Fig. Q3] demonstrates this effect by showing the 
P{s) distribution for the circular QD at locus D in the 
{CCCC} plate, for both the DP and PZ interactions in- 
cluded with (3dp — (ipz = 10. The agreement with 
the GUE statistics is excellent, in contrast to case D of 
Fig. H (we recall that the {CCCC} and {FFFF} cases 
give similar results). Because the AlGaAs alloy is a weak 
piezoelectric material, the DP coupling shows a stronger 
effect in promoting the chaos, whereas the main action 
of the PZ interaction (in the presence of DP) is to break 
the system's overall symmetry. The explanation for such 
a change in the spectrum statistics is left to Section V. 

B. Rectangular 2DEG 

Chaos is also observed in the calculations for a rect- 
angular 2DEG interacting with the suspended phonon 
cavity. In this section we investigate such nanostruc- 
tures following the procedures previously described. The 
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TABLE II: 
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FIG. 8: The nearest neighbor spacing distribution for the 
{CCCC} plate with the circular QD located at D. Both the 
DP and PZ interactions are included with (3dp = Ppz = 
10. The curves correspond to GOE (dashed) and GUE (dot- 
dashed) . 



FIG. 7: Top panel: A 3 (Z) statistics at point A of the {CFFF} 
plate for various radii of the QD: R = 200 nm (filled dot), 
300 nm (star) and 400 nm (open dot); Bottom panel: same 
statistics for various m*/m e ratios: 0.2 (filled dot), 0.4 (star), 
0.6 (open dot), and 0.8 (cross). /3dp=W. 



obtained statistics are summarized in Table II for the in- 
teractions, either DP or PZ, taking into account individ- 
ually. The calculations were made for the same phonon 
cavity considered throughout Section liV Al but now sup- 
porting a square QD of sides 400 nm and thickness equal 
to the circular case. 

Representative results of the P(s) distribution are 
shown in FigEJ which illustrates the chaotic behavior of 
the {CCFF} phonon cavity through cases A to E. Here 
too, the global symmetries of the system depend on the 
combination of the symmetry axes of the plate with the 
position of the square QD. From extensive simulations 
we verified that (see Table II) whenever the full problem 
has only one global symmetry axis, either 8 X , Sy, S x —y 
or S x+y , the resulting spectral statistics corresponds to 
the superposition of two uncorrelated GOEs, contrasting 
with the case of the circular QD-NEMS (see Table I). If 
there are no overall symmetry axes, the statistics is that 
of the GOE type. Finally, the spectrum is regular if there 
are at least two global symmetry axes, namely, the posi- 
tion A for the {CCFF}, {CCCC} or {FFFF} boundary 
conditions. 



Despite the fact that circular and rectangular QD- 
NEMSs display chaotic features, it is important to em- 
phasize that the GUE statistics never occurs for the rect- 
angular 2DEG coupled to a rectangular phonon cavity. 
This is, therefore, an effect that results from the inter- 
play between the cylindrical and rectangular symmetries 
in the circular QD-NEMS. We discuss this phenomenon 
in detail in the next Section. 

The dependence of the spectrum statistics on the geo- 
metrical and material parameters is, nonetheless, similar 
to that observed for the circular QD. Specifically, heavier 
in-plane electron effective masses, lighter and softer ma- 
terials and larger and thinner quantum dots favor the ap- 
pearance of chaos. The {CFCF} plate requires stronger 
interaction strengths than the other boundary conditions 
to give rise to a chaotic spectrum. 

Finally, when the system has one global symmetry axis 
and both the DP and PZ interactions act simultaneously, 
the spectrum statistics changes from two uncorrelated 
GOEs to a single GOE. This effect is illustrated in Fig 
^]for the rectangular 2DEG centered at points A (right 
panel) and D (left panel) of the {CFFF} cavity. A com- 
parison with Table II evidences the aforementioned trans- 
formation. On the other hand, when the nanostructure 
displays either regular or chaotic (GOE) distributions for 
one of the interactions, the inclusion of the other does not 
alter the original statistics, regardless of the strength of 
the interactions. 
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The Hamiltonian block structure due to the 
phonons 



FIG. 9: P(s) for the {CCFF} boundary conditions and the 
rectangular QD centered at positions A, B, C, D, and E. Only 
the DP interaction is considered (/3d p = 10). The dashed 
lines represent the uncorrelated superposition of two GOEs 
and the continuous lines the regular and GOE cases. 




FIG. 10: P(s) for the {CFFF} boundary conditions and both 
interaction potentials acting together {Pdp = Ppz = 10). 
The rectangular QD is located at positions A (right panel) 
and D (left panel). 



V. DISCUSSION 

It has been shown that distinct geometrical configu- 
rations of the QD-NEMS produce different energy-level 
statistics, in most cases typical of chaotic dynamics. To 
understand this effect we investigate the structure of the 
Hamiltonian matrix of our systems and explain the previ- 
ous results in terms of the underlying symmetries of the 
problem. At the end, we explain the anomalous GUE 
statistics in the light of a more general analysis^SiSS. 
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FIG. 11: (a) A block of the interaction potential matrix 
representing fixed electron quantum numbers and all possi- 
ble combinations for the phonon states. The small blocks d, 
shown in (b), are diagonal matrices whose diagonal elements 
du are all equal. The t's forming the main diagonal in (a), de- 
picted in (c), are also block tridiagonal matrices. Noticeably, 
d and t in (a) are different from the corresponding blocks in 
(c). 

From the Eqs. (|18fl and i|26|l . for the deformation and 
piezoelectric potentials, one verifies that the interaction 
mechanism is mediated by one phonon processes. This 
becomes clear by writing the matrix elements in the basis 
m K= («,n)] 



C I \ ^ K " K ' 



x \pn'{ ni 



Sn'^ n' a -l ■ ■ ■ 8 n » n ' N 



+ 1 



(36) 



where C andX QK " K ' denote, respectively, the appropriate 
coupling constant and the overlap of the phonon mode a 
with the electronic eigenfunctions k" and k'. Notice that 
the Kronecker S's allow only a single phonon transition. 

In this representation we have a very particular form 
for the interaction matrix. Consider the block n x n, 
schematically depicted in Fig. ^](a). It corresponds to 
fixed values for the electron quantum numbers k" and k', 
but embraces all possible configurations for the phonon 
states. The structure of the block n x n is such that the 
outermost block spans all possible states for the phonon 
quantum number n\ = 0, 1, . . . , n. Then, the next inter- 
nal block spans the quantum number ri2, followed by the 
inner blocks 713,714 ... njy. Due to the action of the Kro- 
necker S's in l|36l) . the n x n matrix is block tridiagonal. 
Consequently, the small blocks d (refer to Fig. ^2 (b)) 
must be diagonal, since the interactions are mediated by 
one phonon only. On the other hand, the small blocks t 
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are also tridiagonal (Fig. ^] (c)). Such self-similar ar- 
rangement goes on at all block levels % n 2 ... n jy. 



B. Phonon mode parities 

The reflection symmetries of the phonon cavity lead 
to phonon modes of well defined parity. Such properties 
are examined in this section; for guidance refer to Fig. 
03 For instance, when the boundary conditions at Y\ and 
Y2 are equal, i.e., both C or both F, the modes a have 
either a symmetric (+) or an anti-symmetric (— ) parity 
with respect to S x - The same holds for S y regarding 
the edges X x and X 2 . If the boundary conditions are 
opposite at X\ and X2 and also at Y\ and Y2, then one 
of the main diagonals of the plate, S x - y or S x + y , is the 
only symmetry axis. Therefore the phonon modes a have 
well defined + and — parities about it. Finally, the cases 
{CCCC} and {FFFF} have definite parities about all the 
symmetry axes: S x , S y , S x - V and S x+V . 



C. Circular quantum dot 

In the following we analyze the element I a K >> K > , ap- 
pearing in Eq. i|oT)|) . For the circular quantum dot case 
we have, generalizing Eqs. I|21|) and J57J, 



(a) 



y 

C f (F) 



(b) 



C i (F) 



Ian"*' = f(d) / dxdy g K » K >(r) F a (x,y) 

x {cos[(l" -l')Q]+isin[(l" -l')Q]}. (37) 

Here, g K " K /(r) denotes the product of Bessel functions 
coming from Eq. I|16|) . r = y{x~ x o)' 2 + (y ~ Do) 2 with 
{xo,yo) as the center of the QD, 9 is measured from the 
S x axis and f(d) results from a simple integration along 
the z axis. For the deformation potential F® p (x,y) = 
V 2 J2m n X m {x) Y n (y), whereas for the piezoelectric 
interaction F pz (x,y) = ^^J2 m , n A mn X m(x)Y n (y). 
Notice that the Laplacian is a second order operator, 
therefore the function F® p has the same parity as the 
phonon mode a. On the other hand, F pz results from 
first order derivatives, so it has the opposite parity of a. 

A first examination of Eq. l|o7[) reveals that if I" = V 
the sine function vanishes and the matrix element is real. 
For I" 7^ V, it will be complex, real or purely imaginary 
depending on the system's characteristics. For the sake 
of understanding we shall briefly analyze a representative 
case. Let us assume the same boundary conditions for Yi 
and Y2, then consider two scenarios: the circular QD at 
locus A and different boundary conditions for X\ and 
X 2 (Fig. EI a ))j or locus D regardless of X\ and X 2 (Fig 
E2{k)). In both cases only the phonon parity about the 
S x axis will be relevant for the evaluation of (|37|l . In fact, 
the + (— ) parity of F a about S x leads to a matrix element 
that is real (purely imaginary) for the DP interaction and 
purely imaginary (real) for the PZ interaction. That is a 
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FIG. 12: Circular quantum dot positioned at: (a) locus A 
for the {FCCC} (or {CFFF}) plate; and (b) locus D for the 
{X1X2CC} (or {X1X2FF}) plate. Notice that both cases 
have overall symmetry only about the x axis, i.e., S x - 



consequence of the parity of the sine and cosine functions 
about = together with the parity of F a regarding the 
same axis. 

On the basis of the above analysis and the discussion 
of Section IV Al it follows that if there is one or more 
global symmetry axes in the system, then each matrix 
block n x n of Fig. II If a) can be written as A + iB, with 
A and B originating from the cosine and sine parts of 
the integral in Eq. l(o7|) . respectively. Moreover, those 
matrices are real symmetric and mutually disjoint, that 
is, for A rs 7^ (B rs 7^ 0) necessarily M rs = (A rs = 0). 

Therefore, when a single interaction mechanisms is act- 
ing, we have the following scenarios: 

• If the geometrical configuration of the nanostruc- 
ture is such that there is only one global symme- 
try axis (e.g., locus A or D for plate {CFFF}), 
the ensuing partial symmetry break is enough to 
generate chaos. Moreover, the matrix represen- 



tation of the Hamiltonian H 



Hn+ H 



DP(PZ) 



can be written in blocks of fixed Us as o.;//;/ = 
A;«i' + (i) si & a ( 1 "- 1 ') Br;/, where A and B are dis- 
joint, real and symmetric. Thus, H is completely 
characterized by orthogonal matrices, so belonging 
to the GOE universality class. It is straightforward 
to verify that the present reasoning encompasses all 
the cases of a single GOE statistics listed in Table 
I. 

For locus A and {CCFF} boundary conditions, the 
above structure for H is still valid. However, now 
the system has two symmetry axes, leading to new 
restrictions for the matrix elements. In fact, denote 
by a x <Ty (with a = ±) the a mode parities with re- 
spect to S x and Sy. One finds that the integral over 
the cosine (sine) in Eq. I|37fl is different from zero 

only if \l" — 1'\ is even and the a mode is ++ ( ), 

or \l" — l'\ is odd and the a mode is +— ( — h). 
Such selection rules produce two different families 
of eigenvalues for the problem. For {CCFF} (and 
{FFCC}) each distinct family is chaotic, explain- 
ing the occurrence of two superposed uncorrelated 
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Wigner distributions in the P(s) statistics (for an 
explicit example, see the simpler case of a rectan- 
gular QD in Sec. V-D). 

• For locus A and boundary conditions {CCCC} (or 
{FFFF}) there exists one further global symme- 
try, namely, the equivalence of the x and y direc- 
tions. The extra symmetry prevents the emergence 
of chaos. 

• Finally, in the absence of a global symmetry axis 
{e.g., the quantum dot at C for any boundary con- 
dition, or loci B, C or E for {CFFF}) the Hamilto- 
nian matrix does not separate into real and purely 
imaginary disjoint parts. Hence, it is a full complex 
unitary matrix and the chaotic behavior takes place 
with the system belonging to the GUE universality 
class. 

The last case to be considered is the inclusion of both 
interactions in the Hamiltonian. From the previous dis- 
cussion we know that for a given parity of the mode a 
the DP and PZ potentials lead to exactly opposite types 
of matrix elements. Indeed, if the DP matrix element is 
real (pure imaginary), necessarily that corresponding to 
PZ is pure imaginary (real). Therefore, the Hamiltonian 
H = Hq + Hop + Hpz has a complex matrix representa- 
tion that results in a GUE statistics for the energy levels. 
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FIG. 13: Top: schematics of the 32 x 32 interaction matrix for 
locus D in the {CCFF} nanostructure, ordered aspgai 012 as 



with 



1 or 2 and a; = or 1. Bottom: the transformed 



matrix in a block form. The filled dots indicate the nonzero 
elements. 



D. Rectangular quantum dot 

For the rectangular quantum dot the element T aK >> K > 
can be written as 

1 aK "K' = f{d) j dxdy F a (x,y) 

L X Ly 

x sin^V-^f^] S m[g / 7T^f% (38) 

L X Ly 

where x — xq — L x /2, y = yo — L y /2 and (xo, yo) are the 
coordinates of the center of the QD (for guidance refer 
to Section nTB"|) . From Eq. it is evident that the 

matrix elements are always real numbers. Consequently, 
any chaotic behavior must belong to the GOE class and 
the occurrence of the GUE statistics is ruled out for this 
nanostructure. 

With respect to the quantum numbers, the conditions 
for which the above integral is different from zero are 
again entirely dependent on the global symmetries of the 
system. For instance, if the whole nanostructure has S^, 
as a symmetry axis, then X aK n K i is nonzero for the fol- 
lowing combinations: p" + p' = even and a mode parity 
<j x = +, or p + p' = odd and 172; = —. Similar relations 
hold for q" + q' regarding S y . 

The behavior of the spectral statistics generated by 
the rectangular QD-NEMS can be summarized by the 



following representative situations: for the loci B, C or E 
in the {CFFF} plate there are no global symmetry axes 
and we obtain GOE distributions. For loci A and D in 
the {CFFF} plate, there is one overall symmetry axis 
(Sa;) and the resultant statistics is the superposition of 
two uncorrelated GOE distributions. Finally, for locus A 
in the {CCFF} and {FFFF} (or {CCCC}) plates, which 
contain more than one global symmetry axes, no chaotic 
features are observed. One can verify that all cases in 
Table II follow the same trends. 

In order to visualize the occurrence of the 2-GOE 
statistics, consider the case D in the {CCFF} plate. De- 
spite the fact that the phonon modes a have two sym- 
metry axes, Sa, and S y , only the parity about S^ is a 
global symmetry, due to the position of the QD. As- 
sume then 3 phonon modes, such that a basis state 
is written as \p, q; m, n%, 713), with p,q—l or 2 and 
?i Q =0 or 1. In addition, the 3 phonon mode parities 
with respect to S^ are taken to be {+,—,+}. It re- 
sults in the 32 x 32 matrix schematically represented 
in the top of Fig. ^3 where the filled dots indicate 
the nonzero elements. It is possible to transform the 
original matrix in that shown at the bottom of Fig. 
1131 just by rearranging its rows and columns. By la- 
belling the original rows (from left to right) and columns 
(from top to bottom) as 1, 2, . . . , 32, we obtain the first 
nonzero block of the transformed matrix by perform- 
ing the operation 12 3 4 5 6 7 8 91011 1213141516 -> 
1 2 5 27 13 19 10 9 32 31 6 20 14 28 23 24. A similar proce- 
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dure, i.e., operating over the remaining 17... 32 posi- 
tions, leads to the other nonzero block. Here, the el-ph 
interaction generates chaos in each family of eigenval- 
ues, originating from the two independent blocks. Con- 
sequently, the spectrum of the full matrix gives rise to 
the superposition of two uncorrelated GOE distributions. 
The above analysis is valid for any matrix size. 

Finally, if both interactions act together in a system 
with a single global symmetry axis, say S^, their effect 
is to break the selection rules previously described. This 
happens because F a (x,y) has opposite parities for the 
DP and PZ interactions. As a result the Hamiltonian ma- 
trix does not have a block form and a pure GOE statistics 
emerges from the 2-GOE case, as seen in Fig. 1101 



E. Symmetry operator analysis of the anomalous 
GUE statistics 

So far we have examined the structure of the Hamilto- 
nian matrix to explain the chaotic features exhibited by 
our NEMS. Here, we make a link between our results and 
a more general analysis^ to clarify the appearance of the 
anomalous GUE statistics in our time-reversal invariant 
(TPJ) system. 

As already mentioned, the spectral fluctuations of TRI 
chaotic systems typically correspond to the GOE dis- 
tribution. However, Leyvraz et al— have shown that 
there are exceptions to this rule, which can be inter- 
preted even semiclassicallji2S. Suppose a TRI chaotic 
system that has a discrete point symmetry represented 
by the operator S, then [H, S] = [H, T] =0, where 
H is the Hamiltonian and T the time reversal opera- 
tor. More importantly for the effect, assume also that 
S has two invariant subspaces whose representations are 
complex conjugate of each other. We call them {'I'M} 

and {* (_) }, which are solutions of H^n = . 
Since T*l ±) = [*i ±} ]* = *i T) , it may seem that the 
problem is not TRI because each subspace changes into 
the other under T, therefore causing a GUE statistics 
(notice that the Hamiltonian matrix is complex Hermi- 
tian in this basis). However, this is just an artefact of the 
particular structure of the subspaces. Actually, the full 
Hilbert space is TRI, as can be verified after the simple 
basis transformation $W = i( ±1 - 1 )/ 2 [$C+) ± ^rC-)]/^, 

for which T$„ = . Note also that the Kramers 
theorem2L2& imposes E n + ^ — E n ' . Finally, as pointed 
out in Ref. 28], the present phenomenon is rare because 
often there exists an extra operator V (e.g., the parity 
symmetry operator) for which [H, V\ — P", T 3 ] = 0. This 
operator is responsible for combining the two complex 
conjugate representations of S into an irreducible repre- 
sentation that is self-conjugated, therefore producing a 
GOE statistics. 

Prior to our earliest paper—, the only systems known 
to show such behavior were billiards with three-fold but 
no mirror (parity) symmetries, which have been realized 



experimentally in microwave cavities2£L2ii22. They are 
chaotic by construction (due to their particular geom- 
etry) and have their eigenstates composed by complex 
degenerate doublets (of GUE statistics) and real singlets 
(of GOE statistics). It is possible, however, to establish 
a parallel between our circular QD-NEMS and these bil- 
liards. In our case the electron states, Eq. I|16|l . naturally 
provide the necessary complex representation through 
the angular momentum quantum number I. They are 
divided in singlets, for I = 0, and degenerate doublets, 
for I = ±1, ±2, .... Of course, the original electron states 
as well as the phonon states are regular, but the el-ph 
coupling generates chaos. If, nevertheless, the boundary 
conditions and the location of the QD are such to give 
rise to an global symmetry axis, the energy-level statis- 
tics is of the GOE type due to the ensuing definite parity. 
On the other hand, in the absence of an overall symmetry 
axis (e.g., location C for any plate), no V operator exists 
and GUE statistics arises. 

As a last comment, we recall that in our system the 
original electron degeneracies are destroyed by the inter- 
action with the phonons. Nonetheless, the last behave 
as a perturbation for the electronic spectrum, because 
the energies of the electrons are much higher than those 
for an individual phonon. It is important to mention, 
however, that the occurrence of the doublets is not nec- 
essary for the manifestation of the GUE statistics. Actu- 
ally, even when an additional small perturbation breaks 
that degeneracy, the GUE statistics also arises for each 
split family of eigenstates. It has been confirmed experi- 
mentally by the study of imperfect three-fold microwave 
triangular billiards^. 



VI. CONCLUSION 

We have presented a theoretical study of the 
electron-phonon coupling in nanoelectromechanical sys- 
tems (NEMS) comprised of a suspended dielectric plate 
and a quantum dot on its surface. It is shown that a 
quantum chaotic behavior develops as a result of the el- 
ph interaction, for a wide range of geometrical and ma- 
terial parameters of the QD-NEMS. A method is devel- 
oped to treat this novel class of systems. It associates 
the phonons with the vibrational modes of a suspended 
rectangular plate, for clamped and free boundary con- 
ditions. The electrons are confined to a large QD, of 
either circular or rectangular symmetry, and described 
by the free electron gas approximation. The deformation 
potential and piezoelectric interactions are included non- 
perturbativcly in the model, by calculating the cigencn- 
ergies of the NEMS on the basis of the el-ph states. 

By performing standard energy-level statistics we 
demonstrate that the resulting spectral fluctuations are 
very well described by those of the Gaussian Orthogo- 
nal Ensemble (GOE) or the Gaussian Unitary Ensem- 
ble (GUE). It is evidenced that the combination of the 
phonon mode parities together with the position of the 
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QD determine the overall symmetries of the system, 
whieh arc ultimately the responsible for the distinct 
chaotic features observed. Although, quantum chaos is 
commonly obtained in the system, the GUE statistics 
occurs only in the case of a circular QD-NEMS. It rep- 
resents an anomalous phenomenon, since the problem is 
time-reversal invariant. The fundamental reason for this 
effect lies in the structure of the electronic spectrum, 
which is formed by doublets with I = ±1, ±2, .... In the 
absence of any overall geometrical symmetry, the com- 
plex conjugate doublets transform into each other under 
the action of the time reversal operator, thus simulating 
the behavior of a non-TRI system. 

Finally, calculations are under way to include the ef- 



fects of the electron-electron interaction in the model. 
We conjecture that the same chaotic behavior can also 
arise in this case, because the el-el interaction preserves 
the total angular momentum of the electronic system, 
justifying the previous analysis. 
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